Cluster virial expansion for the equation of state of partially ionized hydrogen plasma 
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We study the contribution of electron-atom interaction to the equation of state for partially 
ionized hydrogen plasma using the cluster-virial expansion. For the first time, we use the Beth- 
Uhlenbeck approach to calculate the second virial coefficient for the electron-atom (bound cluster) 
pair from the corresponding scattering phase-shifts and binding energies. Experimental scattering 
cross-sections as well as phase-shifts calculated on the basis of different pseudopotential models are 
used as an input for the Beth-Uhlenbeck formula. By including Pauli blocking and screening in 
the phase-shift calculation, we generalize the cluster-virial expansion in order to cover also near 
solid density plasmas. We present results for the electron-atom contribution to the virial expansion 
and the corresponding equation of state, i.e. pressure, composition, and chemical potential as a 
function of density and temperature. These results are compared with semi-empirical approaches 
to the thermodynamics of partially ionized plasmas. Avoiding any ill-founded input quantities, the 
Beth-Uhlenbeck second virial coefficient for the electron-atom interaction represents a benchmark 
for other, semi-empirical approaches. 

PACS numbers: 52.25.Kn, 52.20.Hv 
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I. INTRODUCTION 

The thermodynamics of dense plasmas, in particular 
hydrogen, has become an important topic. Models of 
planetary interiors and stars depend on the equation of 
state (EOS) of the most abundant elements Reli- 
able EOS data for hydrogen are indispensable for the 
planning and conduction of inertial confinement fusion 
experiments Q. Within the quantum statistical ap- 
proach to the EOS, numerical techniques like density- 
functional theory and quantum molecular dynamics sim- 
ulations have been elaborated Q. Alternatively, ana- 
lytical approaches have been derived using many-body 
perturbation theory [J] that yield, e.g., rigorous results 
for limiting cases, such as the low density limit, that are 
hardly accessible by numerical methods. 

A rigorous result for systems with short-range inter- 
action is the Beth-Uhlenbeck formula that expresses the 
second virial coefficient in terms of the scattering phase- 
shifts and the bound state energies @. Some gener- 
alizations have been performed to obtain results for a 
larger range of densities by including medium effects, see 
0| f° r charged particle systems or Q for nuclear sys- 
tems. Another particularly successful method to inves- 
tigate strongly correlated many-particle systems is the 
cluster-virial expansion Q . In addition to the elementary 
constituents, the different bound states (cluster) are con- 
sidered as new reacting species in thermal equilibrium. 
The thermodynamic properties are expanded in terms of 
the fugacities of the various components of the system 
and cluster-virial coefficients are introduced. In the case 
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of nuclear systems [9Hll(, one can consider the nucleons, 
deuterons (two-body cluster), and alpha-particles (four- 
body cluster) as well as possible further nuclei as the com- 
ponents of the system. For ionic plasmas, electrons, ions, 
atoms, and molecules are such components. Another ex- 
ample is the electron- hole exciton system [l2|, [u| . 

The cluster-virial expansion represents a "chemical" 
picture in the sense that the virial is expanded in orders of 
the fugacities of the different components (single-particle 
states, bound states) in the system. In contrast, the 
traditional virial expansion is a "physical" picture, i.e. 
the fugacities of elementary particles, such as electrons 
and nuclei, are the expansion variables. In the physi- 
cal picture, bound states appear in higher order virial 
coefficients; their treatment involves sophisticated math- 
ematics. On the other hand, bound states are naturally 
included already in the lowest order of the chemical pic- 
ture. For special parameter values, bound state forma- 
tion gives the leading contribution, e.g. in atomic or 
molecular gases. The chemical picture accounts for these 
main terms already in the lowest order of the virial ex- 
pansion, whereas within the physical picture we have to 
consider higher orders of the expansion to identify the 
leading contributions. 

We consider a partially ionized hydrogen plasma, 
which consists of three components as electrons (e), ions 
(i), and hydrogen atoms (a). The formation of heavier 
clusters, such as molecules or molecular ions, can also be 
included but this is not considered in the present work. 
Restricting ourselves to these three components, the rel- 
evant interactions arc the elementary Coulomb interac- 
tion (e — i, e — e, i — i) and the more complex interaction 
(e — a, i — a, a — a) with the atoms. The Coulombic con- 
tributions to the EOS have been investigated extensively 
elsewhere, see Ref. @, and will not be given here. 
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The treatment of electron-atom interaction in EOS 
studies is still an open question. A widely used method 
consists in constructing an effective electron-ion poten- 
tial. The pseudopotential method allows to include 
medium effects and thereby to enlarge the region of appli- 
cability of the cluster-virial expansion. Such pseudopo- 
tentials should reproduce calculated phase-shifts and ex- 
perimental scattering data. The disadvantage of the 
pseudopotential method is that it depends on a lo- 
cal, energy-independent potential that can be introduced 
only in certain approximations to replace the energy de- 
pendent three-particle T-matrix. Well-known examples 
are the Buc king ham potential [HI, see also more recent 
approaches fl5j. Other semiempirical methods, such as 
the excluded volume concept [1 61 ] discussed below, de- 
pend on free parameters that lack a proper quantum- 
statistical foundation. 

In this work, we overcome the shortcomings of the 
pseudopotential method by deriving results for the 
electron-atom virial coefficient that are based on ex- 
perimental data for the electron-atom scattering cross- 
section. Therefore, our results represent benchmarks for 
the electron-atom contribution to the EOS, only limited 
by the accuracy of measured scattering data. We com- 
pare our results to different pseudopotential calculations. 
As an example, we construct a separable pseudopotential 
for the electron-atom interaction, and consider medium 
effects such as self-energy and Pauli blocking. 

We present results for the EOS, in particular the pres- 
sure, composition, and chemical potentials for tempera- 
tures in the range of T < I0 5 K and for densities in the 



order of n < 10 21 cm 3 , where : 



-.total 



-.total 



denotes 



the total number density of electrons or ions, respectively, 
that are found in free or bound states, and obey the neu- 
trality condition. For these parameter values, the hydro- 
gen plasma is moderately coupled with the plasma pa- 
rameter T = e 2 /(47re fc B rd) < 1, where d = (S/inn) 1 / 3 
is the average interparticle distance. Effects of Fermi de- 
generacy are weak for O = k B T/Ep > 1, with the Fermi 
energy E-p — h 2 (3Tr 2 n) 2 ^ 3 /2m c where, more consistently, 
only the density of free electrons should be taken. With 
increasing density and/or decreasing temperatures, the 
system becomes degenerate. 

The present work is organized as follows: In Sec [Til we 
briefly review the cluster-virial expansion and the Bcth- 
Uhlcnbeck formula for the second virial coefficient. Sec- 
tion Hill contains the calculation of the scattering phase- 
shifts for the electron-atom system, both via experimen- 
tal differential cross sections and phase-shifts from appro- 
priate pseudopotentials. In Sec. IIII El the phase-shifts 
are used to calculate the corresponding second virial 
coefficient. Based on these results, corrections for the 
pressure and the chemical potentials are calculated in 
Sec. IIV1 Sec. [V] contains the presentation of the general- 
ized Beth-Uhlenbeck approach and the calculation of the 
density dependent second virial coefficient. Conclusions 
are drawn in Sec. I VII 



II. CLUSTER VIRIAL EXPANSION AND THE 
BETH-UHLENBECK FORMULA 

The canonical partition function of an interacting 
many-particle system at temperature T, volume V, com- 
posed of N c particles per species c (c = e, i, a) carrying 
spin s c , reads 



Z can (T,V,N c ) = Tr{exp(-/3ff)} 



(1) 



with 3 = l/k B T. Here and in the following, the spin 
of particles of species c is implicitely taken as s e = 1/2 
for electrons, Si = 1/2 for protons and for atoms in the 
singlet state (antiparallel spins of electron and proton) 
Sa.singict = 0, and s Q:tri pi ct = 1 for the triplet state (par- 
allel spins of e and i). We neglect hyperfine splitting of 
the atomic levels. 
The Hamiltonian 



c j = l 



F 



2m 

N c N, 



cd j — 1 k=l 

contains, besides the kinetic energy of each particle, the 
mutual interaction, represented by the two-particle in- 
teraction potential V cc i(fj — f%)- The prime indicates 
that self-interaction is excluded, and the energies for each 
component are gauged to E^ if the particles are at rest. 
From a relativistic approach, Ec is given by the rest 
mass containing the binding energy. In our nonrelativis- 
tic approach, we choose the gauge relatively to the rest 
mass of the elementary particles so that Ec is the bind- 
ing energy of the composites. 

Using the Mayer cluster expansion [ljj , we arrive at the 
virial expansion for the free energy F — -^T lnZ can 
valid for short-range potentials, 

F{T, V, N c ) = F id (T, V, N c ) - k B Tv{ ^ n c n d b cd 

cd 

+ ^ n c n d n e b c de + ...}. (3) 



cde 



In 



AgAf c 
Vgc 



1 + BE C ' 



(0) 



is 
non- 



{2ith 2 /k B Tm c ) 1/2 de- 



F ld (T,V,N c ) = k B TJ2 c N t 
the free energy for the classical ideal (i.c 
interacting) system, where A 

notes the thermal wavelength of species c, and n 
N c /V the particle number density of the component c, 
g c = 2s c + 1 is the spin degeneracy factor. The expan- 
sion coefficients b c d and b cc i e are the second and third 
virial coefficients, respectively. They are determined by 
the interaction, but also by degeneracy terms. 

Having the virial coefficients at our disposal, we can 
easily derive the thermodynamic properties of the system 
under consideration. E.g. for the pressure and the chem- 
ical potential the following expressions are found using 
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the standard thermodynamic relations: 

p(T, V, N c ) = p id (T, V, N c ) - fc B r{ J2 n cn d b cd 

cd 

+2^2n c n d n e b cde + ...X , (4) 

cde 

Hc(T, V, JV C ) = Mc,id(T, V, N c ) - k B T[2 £ n d 6 cd 

d 

+3^2n d n e b cde + ...} , (5) 

where p id (T, V, iV c ) = k B Tj2 c ^ and (j, Cl a(T, V, N c ) = 
ksTln { ^Yg j arc ^ e ideal P ar ts of the pressure 

and the chemical potential, respectively (for the hydrogen 
atom the degeneracy factor is g a = 4). Note that in 
relativistic approaches the chemical potentials are gauged 
including the rest mass of the constituents, as discussed 
for the Hamiltonian ||3J|. 

Working with the grand canonical ensemble, it is con- 
venient to introduce the fugacities z c = e^ c ~ Ec ) anc j 
to expand the grand canonical potential with respect to 
the fugacities [13], see also App. [5] 



-pV 



-k B TV 
^ b cd z c z d 



(0) 



where b cd = b cd g d / A d are the dimensionless second virial 
coefficients. In the low density case n c A;? <C 1, the rela- 
tion z Cj id = n c A^./g c results for the ideal system. Since 
the different ensembles are equivalent in the thermody- 
namic limit, the fugacity expansion is shown to be equiv- 
alent to the expansion with respect to the particle den- 
sities after substitution of the variables; for details see 

Chemical equilibrium for a reaction v a A + v B ±=t 
v c C between the components A, B, C gives the relation 
^afJ-a + ^bfJ-b = v c (i c . This way, the thermodynamic po- 
tentials finally depend only on the total densities of the 
constituents, or their chemical potentials, since the to- 
tal number of the constituent particles is conserved. The 
densities of the composites or their chemical potentials 
can be eliminated, using a mass action law or a Saha 
equation. 

Note that it is possible to derive the cluster virial ex- 
pansion in a systematic way, starting from a quantum 
statistical approach The spectral function of the 

elementary particle propagators is related to the self- 
energy. Within a cluster decomposition of the self-energy, 



the contribution of the different clusters can be identified 
considering partial sums of ladder diagrams. The first- 
principle approach gives a consistent introduction of the 
chemical picture avoiding double countings, and allows 
for systematic improvements. 

In this paper we are concerned with the evaluation 
of the second virial coefficient that describes the non- 
ideality corrections in lowest order with respect to the 
densities. According to Beth and Uhlenbeck 0, H, [HI] , for 
a central symmetric interaction potential the following 
formula has been derived, which relates the second virial 
coeffient to the energy eigenvalues of the two-particle 
bound state \n£) and the scattering phase-shifts rjf d (E) 
describing the two-particle scattering states, 



9d 



5 cd b^>+bT a(i + bf d 



where 



r(o) _ I 2 5 / 2 (ideal Bose gas) 
cd ~ i -2" 5 / 2 (ideal Fermi gas) 



(7) 



(8) 



is the second virial coefficient for the ideal quantum gas. 

In this work we calculate the second virial coefficient 
for the e — a contribution. The spin degrees of freedom 
of the bound electron and the scattering electron give 
rise to a singlet (antiparallel electron spins) and a triplet 
(parallel electron spins) scattering state; the proton spin 
contributes a spin degeneracy factor gi — 2. The total 
second virial coefficient b ae — b ea is defined with the total 
density of atoms and electrons. It is decomposed into the 
singlet contribution and the triplet contribution, so that 
b ae = b s ™ glct + &Jj" plct . For convenience, we introduce the 
dimensionless coefficient b that appears in the fugacity 
expansion as b ea = j^-b ea , b ae = j^b ae . Note that b ea is 

no longer symmetric with respect to a change of indices, 
instead, we find b ea = 2 ^ b ae . Since m e <C m a 

we have Ag a /A^ w 1, from which follows that the di- 
mensionless second virial coefficient is again the sum of 
the corresponding dimensionless singlet and triplet coef- 
ficients, b ae = b s ™ elct + ^ r e iplct . 

The bound part of the second virial coefficient for the 
singlet state b^° nnd ' singlct has the following form: 



Abound, singlet 



1 oo 



n>e+i 



(9) 



where £ denotes the angular momentum of the two- 
particle system. The scattering part of the second virial 
coefficient b s JL consist of the singlet and triplet parts: 



4 



^sc, singlet 



-/3E 



dE 



(10) 



^.triplet = i ^(2£. 



1) 



e -£E " ae.triplet (s) ^ 



(11) 



r 



At this point, we would like to make a short note re- 
garding different forms of the Beth-Uhlenbeck formula 
Eq.© that can be found in the literature. We use for- 
mula Eq.0, which has been derived orig inally by Beth 
and Uhlenbeck || , see also Refs [13, fig . After partial 
integration of the Eq. (J7J) one obtains another form for 
the Beth-Uhlenbeck formula, see e.g. [1, Q, where the 
scattering phase shift arises instead of its derivative, and 
from the integration an additional term —rjg (0)/iT ap- 
ears, which is sometimes condensed into the bound part 



Coming back to the partially ionized plasma, the virial 
expansion is diverging for the long-range Coulomb inter- 
action. This refers to the e — e, e — i, i — i contributions. 
Partial summations lead to convergent results, and the 
expansion of the thermodynamic potentials contains also 

1/2 r — i 

terms with n c and logn c , see [6(. The contribution of 
the scattering and bound parts of the second virial co- 
efficnt for atom-atom interaction was calculated in Ref. 
[l9| . In this work, we evaluate the second virial coeffi- 
cient Eq. ([7]) for the electron-atom interaction. 



III. ELASTIC ELECTRON-ATOM 
SCATTERING AND PHASE-SHIFTS 



Experimental data and first-principle 
calculations 



The Beth-Uhlenbeck formula relates the second virial 
coefficient to few-body properties. For the electron-atom 
contribution, the relevant quantities are the phase-shifts 
for the elastic electron-atom scattering as well as the 
possible bound state energies. No direct measurements 
of the electron-atom scattering phase-shift are available 
in the literature, only scattering cross-sections (i.e. the 
modulus of the phase-shift) have been measured. Accu- 
rate data for the angular resolved scattering cross-section 
were obtained by Williams et al. [20j for electron energies 
between 0.58 eV and 8.7 eV and by Gilbody et al. [H[ for 
3.4 eV, see the review in Ref. (22|. The data are shown 
in Figs. [T]and[21 A bound state H~ is measured at en- 
ergy (0.754 ± 0.002) eV (23]. In this measurement, the 
electron affinity was obtained via the threshold en ergy 
for photodissociative formation of ion pairs from H2 [231 ] . 
The ion pair threshold was combined with the ionization 



potential of the hydrogen atom and the bound dissocia- 
tion energy of H2 to obtain a lower bound to the electron 
affinity. This results is in agreement with the theoreti- 
cal value of 0.75421 eV reported by Pekeris [24| for the 
singlet bound state of H~. 

Theoretical calculations for the e — a scattering phase- 
shifts are abundant. The spins of the two electrons are 
combined to a singlet or triplet state, whereas the or- 
bitals are determined by a three-body Schrodinger equa- 
tion and the symmetry condition for fermions. Fre- 
quently used methods are the close-coupling approxima- 
tion (23] , the R- matrix method [26| .direct numerical so- 
lution of the Schrodinger equation [27j , the wave expan- 
sion method [28|,|29|, and variational calculations. Using 
the variational method, Schwartz et al. [3(| obtained the 
phase-shifts for the s-waves (orbital momentum £ = of 
the e — a system) in the singlet and triplet channels. For 
higher orbital moments £=1,2 calculations have been 
performed by Armstaed [3l| and Register 32 1. Compar- 
ing the experimental data with results of numerous theo- 
retical approaches, see Refs [2(| IH, [33J , it was concluded 
that the variational approach is the most reliable and 
most accurate method in reproducing the experimental 
scattering cross-sections. 

In Figs. [T]and Owe show the differential cross-section 
as calculated from the phase-shifts given in Refs [3(J [HJ 
compared to the experimental data from Ref. [2(| ■ Good 
agreement between theory and experiment is observed for 
scattering angles larger than 8 ~ 25°. Deviations below 
25° are due to the neglect of higher orbital momenta 
I > 2 in the variational calculations 201. 



B. Polarization potential models 

Before presenting results for the second e — a virial 
coefficient using the data given above, we discuss the in- 
troduction of pseudopotentials. One of the advantages of 
introducing pseudopotentials is the possibility to describe 
medium effects such as screening and quantum degener- 
acy. Polarization effects at long distances and exchange 
effects at short distances play a key role in the e — a inter- 
action in a dense plasma. Below we use polarization po- 
tential models, which include these many-body features 
to calculate the scattering phase-shifts. Of course, the in- 
troduction of a local and instant pseudopotential for the 
e — a interaction is only an approximation. A rigorous 
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— phase-shift data from Schwartz and Register 
• experimental data 




50 100 150 

scattering angle 9 [deg] 



200 



FIG. 1: Differential electron-atom (e — a) scattering cross- 
section as a function of the scattering angle for electron 



wavenumber k — 0.3 a B , k — 0.4 a B , k 



0.5 a B 1 ; solid 



lines - phase-shift data from the variational method [3(1 1 32 
circles - experimental data [20T |. 



phase-shift data from Schwartz and Register 
experimental data 




50 100 150 

scattering angle 9 [deg] 



200 



FIG. 2: Differential electron-atom (e — a) scattering cross- 
section as a function of the scattering angle for electron 
wavenumber k = 0.6 a B ,k = 0.7 a^ 1 , k = 0.8 a B . 



treatment involves the three-particle T matrix, which is 
non-local due to exchange effects and depends on energy. 
We do not describe this problem in detail here. 

We calculate the phase-shifts on the basis of the wave 
expansion method [29(. This approach yields the phase- 
shift as the solution to the so-called Calogero [28| equa- 
tion, 

d <^ r) = -±U(r)[cosr,f (k,r)J e (k,r) 



— smrjg d (k, r)n e (k, r) 



(12) 



Here, JAk,r) and ne(k,r) are the Rikkati-Bessel func- 
tions [34], U(r) = (2m r / h 2 )V c d(r) , V c d{r) is the in- 
teraction potential and m r is the reduced mass of 
the two-particle system c — d. As initial condition 
we apply ?yf rf (fc,0) = 0, thereby fixing the arbitrary 
phase-offset. The phase-shifts are defined by f]g d (k) — 
linv-^oo rjf d (k, r) and depend only on the wavenumber k. 
As its main advantage, the wave expansion method allows 
to study directly the influence of the interaction potential 
on the phase-shifts. Secondly, equation (|12|) is easier to 
solve than the full Schrodinger equation applying, e.g., 
variational methods. 

As the electron-atom interaction potential we consider 
the Buckingham model [3] 



V ea (r) = - 



2(r 2 +r 2 ) 2 



(13) 



a = 4.5 <2g is the hydrogen atom polarizability (here and 
henceforth we measure distances in units of the Bohr ra- 
dius, a B = 0.529 A). The cutoff radius r Q is used to 
regularize the behaviour at small distances. Its value 
given in the literature [35[ is r = 1.4565 as- However, 
we suggest here the use of a different value ro = 1.033 ae, 
which yields the correct H~ ion ground state energy 
Eq = — 0.7542 eV as the eigenvalue of the effective ra- 
dial Schrodinger equation [6j. At large distances r ^> tq 
the Buckingham potential describes the typical 1 /r 4 be- 
haviour of the electron potential energy in the field of the 
polarizable H-atom. 

Secondly, we employ the polarization potential model 
that was suggested for semiclassical plasmas in Ref . [l5[ : 



V ea {r) 



where A 2 



2r 4 (l~4\ 2 e /rl) 



-Br 



(1 



-Ar 



Br) 

2 



(1 + Ar)) , (14) 



which is a first order non-linear differential equation. 



1a!7^)/(2a 2 ), b 2 = (i - 

y/l - 4X 2 e /r 2 D )/(2Xl), and A e = A c /2tt is the electron 
thermal de Broglie wave-length, m — [eo k^T /n^e 2 ] 1 / 2 
is the Debye radius due to electrons. This model de- 
pends on two parameters A e and ro- At large distances 
the RDO potential is weaker than the Buckingham model 
(fT"3|) due to screening. The strength of the model at short 
distances is given by — ae 2 /8\ 4 . . Fixing rr> — 4.84 as and 
A e = 0.62 ae, the H~ ground state energy is found at the 
correct energy. 

At this point we want to remark that the bound state 
occurs only in the singlet scattering channel, due to 
the strongly repulsive exchange interaction in the triplet 
channel. However, neither the Buckingham model, 
Eq. fT5)l. nor the RDO model, Eq. (HU), take into ac- 
count this exchange effect. A convenient method to over- 
come this problem consists in using a separable potential 
36, 37]. This will be discussed in the following subsection 

ITntl 
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C. Separable potential method 

Separable potentials have been used extensively in nu- 
clear physics to parametrize the nucleon-nucleon interac- 
tion [36fl . It can be shown that any interaction potential 
can be approximated by a sum of separable potentials 

We characterize different channels by spin and angular 
momentum and assume that there is no coupling between 
different channels. We consider a rank-two separable po- 
tential in momentum representation to describe attrac- 
tion at long distances and repulsion at short distances, 

V(p,p') = Xiwi{p)wi(p') + \ 2 w 2 (p)w 2 (p') (15) 

where wi(p), w 2 (j)) are Gaussian form factors Wi(p) = 
cxp (— p 2 /bf), and Ai, bi are the strength and interac- 
tion range, respectively. We determine these parameters 
by fitting the phase-shifts obtained from the separable 
potential to the experimental data. Thereby we exploit 
the definition of the phase-shifts as the argument of the 
T-matrix, 



TABLE I: Parameters of the rank-two and rank-one separable 
potentials (Ai, A2, in Hartree, and 61, 62 in as), the scattering 
length (a), the effective range (R), the binding energy (Eo) 
for the singlet (S = 0) and the triplet (S = 1) channels of 



Ai A2 



61 



b'2 a/a,B R/cib Eq, Hartree 



rank-two separable potential 



s = 
s = 1 



-25.4 10 0.8 0.787 5.965 
37 40 0.8 0.787 1.97 



-0.06899 



rank-one separable potential 



S 

5: 



-45.4 0.4705 
77.67 0.9 



5.965 3.32 
1.77 1.1 



-0.0474 



the bound state properties, can be reproduced in certain 
approximation by separable potentials. We expect that 
increasing the rank of the potential, the experimental 
values for the two-particle properties are better realized. 



tan?y(p) = 



(16) 



D. Phase shifts data 



The calculation of the T-matrix for the separable po- 
tential, which amounts to the solution of the effective 
Schrodinger equation 



P 

2m. 



(17) 



p' 



is detailed in the App. [Bj 

Fig. [3]shows the phase-shifts for the singlet and triplet 
scattering channels obtained by the separable potential 
in comparison with data of Ref. [30j j . The best fit pa- 
rameters are summarized in Tab. HI Column 4, 5 and 
6 give the scattering length a and the effective range R 
from the effective radius theory, and the binding energy 
Eq used to fix two of the parameters Ai and bi, the re- 
maining two being fixed directly by comparison to the 
experimental phase-shifts. The effective radius theory 
was applied e.g. in Ref. [3^| to calculate the influence 
of atomic and molecular contributions to the EOS of hy- 
drogen plasma. However, the effective radius theory is 
limited to s-wave scattering, and therefore to low ener- 
gies. Furthermore, in Ref. [39 ]. the spin dependence of 
the e — a interaction is neglected. 

Choosing X 2 = and b 2 — in Eq. (fTS"]) . a rank-one 
separable potential is obtained. Parameters of a rank- 
one separable potential are given in the Tab. [J (bottom 
part). Using the properties of the T-matrix (see Refs. 
Is big])) we found the binding energy E of the H~ ion for 
both versions of the separable potentials, given in the last 
column of Tab. |TJ The experimental value of the binding 
energy is E Q = -0.0277 Hartree (-0.7542 eV) [23]. The 
two-particle properties, i.e. scattering phase-shifts and 



Schwartz data 
RDO model 

RDO model (bound state) 
Buckingham model (r () -l .4565a B ) 
Buckingham model (r -1.033a B ) 
rank-one separable potential 
rank -two separable potential 

triplet " 




wave number k [a B ] 



FIG. 3: Electron-atom (e-H) scattering phase-shifts as a func- 
tion of wave number k for I = 0. Shown are experimental data 
of Schwartz 30] for the singlet and triplet channels, model 
calculations using the Buckingham [Eq. H13p ] and the RDO 
pseudopotential [Eq. H14p ] at different parameters, as well as 
different separable potentials. 



Using the potentials (fT3"]) and (fT4"| , the Calogero equa- 
tion (fTl?]) is solved. The s-wave scattering phase-shifts 
obtained in this way are plotted as a function of the 
wavenumber k in Fig. [3] We compare our calculations 
to the experimentally validated data by Schwartz et at, 
employing different choices for the cutoff parameter tq in 
the case of the Buckingham potential, and ro and A e for 
the RDO potential. 
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FIG. 4: Electron-atom (e-H) scattering phase-shifts as a func- 
tion of wave number k for 1 = 1. 
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FIG. 5: Electron-atom (e-H) scattering phase-shifts as a func- 
tion of wave number k for i — 2. 

At k = 0, the singlet phase-shifts from [3(| tend to 
r/o = 7r, corresponding to one bound state as follows from 
the classical Levinson theorem ^(j> *?(0) = n7T ( n i s num- 
ber of bound states). The polarization potential method 
gives a bound state if the screening parameters that re- 
produce the correct H~ binding energy are applied (red 
dashed curve for Buckingham potential, ro = 1.033 ae 
and blue square curve for RDO potential, rjj = 4.84 ob 
and A e = 0.62 ae). Taking the original screening param- 
eters in both models, we find vanishing phase-shifts at 
k = 0, i.e. no bound state. The solid line and dash- 
dotted line correspond to phase-shifts for rank-two and 
rank-one separable potential, respectively. The phase- 
shifts from the rank-two separable potential fully coin- 
cide with the experiments, whereas the results for the 
rank-one separable potential deviate at high values of k 
(respectively E). 



We consider the s-wave scattering phase-shifts in the 
triplet channel in Fig. [3J At zero-energy the phase-shift 
starts off at ir. Having in mind that the effective interac- 
tion between the electron and atom in the triplet channel 
is repulsive, this behavior obviously contradicts Levin- 
son's theorem that predicts the scattering phase-shift to 
increase by n with every occurring bound state. To re- 
solve this inconsistency, the classical Levinson theorem 
for one-body problems has been generalized for the scat- 
tering on a compound target by Rosenberg and Spruch 
[4ll ] . The generalized Levinson theorem states that the 
phase-shift at vanishing energy is 77(0) = (rtp au ii + n)ir, 
where np au ii is the number of states from which the par- 
ticle is excluded by the Pauli principle. np au ii is de- 
fined by the number of nodes in the one-particle wave 
function. Application of the generalized Levinson the- 
orem to the electron-hydrogen triplet scattering shows 
the one-particle wave function has one node, that means 
the zero-energy triplet phase shift is a nonzero multiple 
of 7T [13] • Since the triplet electron- hydrogen wave func- 
tion is spatially antisymmetric, the equivalent one-body 
wave function is orthogonal to the hydrogenic ground- 
state function and must have at least one node. Thus, 
our result for the behavior at zero-energy of the triplet 
phase-shift does not predict a triplet bound state H~ . 
This agrees with previous investigations of scattering of 
electrons on hydrogen atoms at low energies, where the 
scattering length was defined under the assumption that 
the negative hydro gen ion H~ can be formed only in the 
singlet channel [43r - l47j ] . In our calculations of the second 
virial coefficient we consider only the singlet bound state 
of if". 

Figs. H] and E] show the calculated phase-shifts for 
p and d waves, respectively. For comparison, the data 
from variational calculations [12] are also presented. The 
phase-shift is zero (no bound states) at the origin and 
increases monotonically. The phase-shifts for I = 1 and 
i = 2 are very small in comparison to the s-wave data 
for the low-energy range. In general, terms from higher 
orbital moments are negligible for low-energies. Hence, 
we only apply phase-shift data for I = 0, 1, 2 in further 
calculations. 



E. Results for the second virial coefficient 

The good agreement between experimental cross- 
sections and the variational scattering phase-shifts from 
Refs. [13, [HI allows us to use the latter as "experimen- 
tally confirmed" data for calculations of the second virial 
coefficient using the Beth-Uhlenbeck formula (J7J). The 
phase-shifts obtained from pseudopotential models (fl"3"|) 
and Q14[) are not applied here to calculate b ae . 

Tabs. |TT] and IIIII show results for the normalized sec- 
ond virial coefficients ^ c / inglet (T), fobovmcMinglet^ for 
the singlet channel and b^ tTlpiet (T) for the triplet chan- 
nel, respectively. The second, third and fourth columns 
of both tables present data for the contribution of s, p 
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and d waves to the scattering part of the second virial 
coefficient for singlet and triplet channels, respectively. 
Higher order contributions are small and negligible for 
this temperature range. The singlet bound part of the 
second virial coefficient for the singlet channel is shown 
as the fifth column of Tab. Hl| T he binding energy is 
taken as E e B H = -0.7542 eV~||. The full scattering 
part of the second virial coefficient is shown in the sixth 
column of both tables. The last column of the Tab. IIIII 
presents the results for the full second virial coefficient: 

b ae (T) = ^ c / in S lot (T) + 6^ triplet (T) + gbound,singlet( T ) < 

We find that the scattering contributions to the second 
virial coefficient increases with temperature in contrast 
to the bound state contribution. In Tabs. ITT1 and ITTT1 the 
s-wave contribution to the second virial coefficient is the 
dominant term, p-wave and d-wave contributions are of 
the order of few percent. 

IV. EQUATION OF STATE AND 
THERMODYNAMICS 

A. Composition 

The virial expansion allows to determine a thermody- 
namic potential that gives all thermodynamic variables. 
We discussed the free energy F(T, V, N c ) or the grand 
potential —pV — J(T,V,fi c ). Because of reactions in 
the system, the particle numbers of the different compo- 
nents are related by the chemical equilibrium conditions 
so that the number of independent variables is reduced. 
In the case of a hydrogen plasma considered here, we start 
from the particle numbers of free electrons, free ions, and 
atoms, disregarding heavier clusters. The atomic density 
is related to the free electron and ion density due to the 
Saha equation that follows from the equilibrium condi- 
tion fj, e + fj,i = [i a . The remaining two particle num- 
bers, 7V* ot = N c + N a with c — e,i, will coincide if a 
charge-neutral plasma is considered so that we end up 
with only one particle number TV = Nl ot = iV* ot for a 
charge-neutral hydrogen plasma in chemical equilibrium. 

To derive the composition from the chemical equi- 
librium condition we express the chemical potentials in 
terms of the densities, see Eq. ([5]). In lowest order of the 
cluster virial expansion, the ideal Saha equation 

^=nl° t Ale X p[pI efi (n e ,T)] , (18) 

is obtained for the degree of ionization a — n e /n t ° t 1 
where I°£(n e ,T) = \E° a \. 

We will not discuss here the more general expressions 
where the excited states and higher clusters are included 
Q. The thermal wavelength for the atom was approxi- 
mated by the thermal wavelength for the ion. 

Taking the non-ideal terms into account, e.g. accord- 
ing to a virial expansion, the composition follows from a 
Saha equation with shifted energies |4J 

/^(nfSnfST) = \E° a \ - A a + A e + A, . (19) 



The energy shifts A c of the different components can be 
obtained from density expansions. As an approximation 
we take the Debye shift A e = Aj = — ne 2 /2 due to the 
Coulomb interaction between the charged particles, k = 
[rig ree e 2 /eo fceF] 1 / 2 is the inverse Debye screening length. 

1/2 

These terms are of the order n e ■ The bound energy 
shift A a is not taken into account here because it is of 
higher order in density. 

In Fig. [S]we plot the solution of the Saha equation (fig)) 
in dependence of the total electron density for temper- 
atures T = 15000 K, 20000 K, and 30000 K. The degree 
of ionization is decreasing with increasing density due to 
formation of bound states. The effective bound state ion- 
ization energy / is lowering due to plasma screening. 
Ultimately, this leads to the Mott effect, i.e. the non- 
thermal ionization at high densities, due to the lowering 
of the ionization threshold, leads to the abrupt increase 
of the ionization degree, see also Refs. [H, @|. We refrain 
from giving an exhaustive description of the Mott effect 
including more sophisticated analysis of the shifts and 
restrict ourselves only to the general behavior of the ion- 
ization degree. Note, that the virial expansion can only 
be applied to the low density range where the corrections 
are small. 




total electron density n [cm ] 



FIG. 6: Degree of ionization as a function of the total electron 
density for three different temperatures, T — 1.5, 2, 3 • 10 4 K. 



B. Chemical potential, free energy, pressure 

To discuss the contribution of electron-atom scatter- 
ing to the chemical potential ([5]), free energy ((3|), and 
pressure (0]) we rewrite the definitions as 

Afi ea = fi e {T, nc) - fif d = -2k B Tn a b ea (T) , (20) 
AF ea = F(T, n c ) - Fid = p(T, n c )V - Pid V 

= -2k B TVn e n a b ea {T) . (21) 

The remaining contributions to the second virial coeffi- 
cient due to the other combinations of components will 
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TABLE II: The singlet scattering and bound parts of the second virial coefficient for e — H interaction. Contribution of different 
partial waves and bound state contribution are given. 



T, K 


6ae' smget ! s-wave 


6ae' Blng et , p-wave 


ba/ inslet , d-wave 


£bound,singlet 


b s ^ slet full 


5000 


-0.0499 


0.0012 


0.0007 


1.4401 


1.3922 


6000 


-0.0541 


0.0015 


0.0008 


1.0756 


1.0239 


7000 


-0.0578 


0.0017 


0.0010 


0.8732 


0.8181 


8000 


-0.0611 


0.0018 


0.0012 


0.7468 


0.6887 


9000 


-0.0642 


0.0020 


0.0013 


0.6613 


0.6005 


fOOOO 


-0.0670 


0.0021 


0.0015 


0.6000 


0.5366 


flOOO 


-0.0696 


0.0022 


0.0016 


0.5541 


0.4883 


f2000 


-0.0720 


0.0022 


0.0018 


0.5185 


0.4505 


f3000 


-0.0743 


0.0023 


0.0019 


0.4902 


0.4202 


14000 


-0.0765 


0.0023 


0.0021 


0.4672 


0.3952 


15000 


-0.0785 


0.0024 


0.0022 


0.4481 


0.3742 


20000 


-0.0873 


0.0024 


0.0029 


0.3873 


0.3054 


30000 


-0.1004 


0.0022 


0.0043 


0.3347 


0.2409 


40000 


-0.1099 


0.0019 


0.0057 


0.3111 


0.2089 


50000 


-0.1172 


0.0016 


0.0072 


0.2978 


0.1894 


60000 


-0.1230 


0.0015 


0.0088 


0.2892 


0.1765 


70000 


-0.1276 


0.0015 


0.0105 


0.2833 


0.1678 


80000 


-0.1312 


0.0017 


0.0125 


0.2789 


0.1620 


90000 


-0.1340 


0.0022 


0.0148 


0.2755 


0.1584 


100000 


-0.1361 


0.0028 


0.0172 


0.2728 


0.1568 



TABLE III: The triplet scattering part of the second virial coefficient for e — H interaction. Contribution of different partial 
waves, last coloumn: full second virial coefficient b ae . 



T, K 


6 s a c e ' triplct , s-wave 


KT iplct , P-wave 


^.triplet ^ d _ wave 


5* r e iplct full 


bae full 


5000 


-0.0548 


0.0121 


0.0024 


-0.0402 


1.3519 


6000 


-0.0604 


0.0147 


0.0029 


-0.0426 


0.9812 


7000 


-0.0654 


0.0174 


0.0033 


-0.0446 


0.7735 


8000 


-0.0702 


0.0200 


0.0038 


-0.0462 


0.6425 


9000 


-0.0746 


0.0227 


0.0043 


-0.0475 


0.5529 


10000 


-0.0788 


0.0254 


0.0048 


-0.0485 


0.4880 


11000 


-0.0828 


0.0281 


0.0053 


-0.0494 


0.4389 


12000 


-0.0866 


0.0307 


0.0057 


-0.0501 


0.4004 


13000 


-0.0902 


0.0333 


0.0062 


-0.0506 


0.3695 


14000 


-0.0937 


0.0360 


0.0067 


-0.0510 


0.3442 


15000 


-0.0970 


0.0386 


0.0071 


-0.0512 


0.3230 


20000 


-0.1121 


0.0512 


0.0094 


-0.0514 


0.2539 


30000 


-0.1367 


0.0742 


0.0139 


-0.0485 


0.1924 


40000 


-0.1567 


0.0943 


0.0183 


-0.0440 


0.1648 


50000 


-0.1735 


0.1117 


0.0224 


-0.0393 


0.1501 


60000 


-0.1882 


0.1270 


0.0264 


-0.0347 


0.1418 


70000 


-0.2012 


0.1405 


0.0302 


-0.0304 


0.1373 


80000 


-0.2128 


0.1525 


0.0338 


-0.0264 


0.1355 


90000 


-0.2233 


0.1635 


0.0372 


-0.0226 


0.1358 


100000 


-0.2328 


0.1735 


0.0403 


-0.0189 


0.1379 
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not be discussed here, see 0, Q. 

As mentioned before, the second virial coefficient can 
be decomposed into the singlet and triplet channel and 
it is given as the sum of scattering and bound state con- 
tributions, 



energy reads 



b ea (T) = b ae {T) 



A 3 
2 



+^° und ' singlct (T)l . (22) 



The various terms are given in Tabs [TT] and IIII1 

In Fig. [7J we plot the e— H scattering contribution 
to the chemical potential Afj, ea = —2kBTn a b ea (T) 
as a function of the total electron number density 
for T = 10000 K and T = 15000 K. In a similar 
way, we treat the e— H contribution to the free en- 
ergy AF ea — —2kBTVn e n a b ea (T) and to the pressure 
Apea = -2kBTn e n a b ea {T) . 



0.5 



0.4 



0.3 



it 

< 0.2 



0.1 

0' 



i experimental phase shifts 
RDO model (bound state) 
Buckingham model (r =1.033 a B ) 



T= 10000 K 




total electron density n [cm 



FIG. 7: Contribution of electron-hydrogen interaction to the 
chemical potential as a function of the total electron density; 
squares represent calculations based on experimental phase- 
shift data, the solid line uses the Buckingham pseudopoten- 
tial, the dot-dashed line uses the RDO model. 



C. Comparison of the results with the excluded 
volume concept 

An alternative approach to evaluate the non-ideality 
term due to the neutral atoms is the excluded volume 
concept [HI]. The excluded volume is defined by the fill- 
ing parameter 77 = ^Trr^n a as the volume that is occupied 
by atoms such that the effective volume available for the 
moving particles is V* = V(l — rj). The atom radius r a is 
an empirical parameter of the order of ob that has been 
fixed in different ways (see Ref. [HI]). Within the ex- 
cluded volume concept, the non-ideality part of the free 



AF e e a x = Fid {T, V*,N e , N u N a ) - F id (T, V, N e , N it N a ) 

(23) 



4 

= -7Trlk B TN a [n e + n l + n a }. 



The corresponding second virial coefficient for the 
electron-atom pair results as 



7TT 

3 



(24) 



It is instructive to note that this expression is equal 
to the classical second virial coefficient within the Bcth- 
Uhlcnbeck approach using a hard-sphere electron-ion po- 
tential with the hard-sphere radius equal to the atom ra- 
dius of the excluded volume concept r a [lj| • It does not 
depend on the temperature of the system. For a typi- 



of 



= 1.0 ae we find U 



■lass 



cal atom radius 
-3.1 x 10" 25 cm 3 . 

In Fig. [5] we show the second virial coefficient for the 
triplet state, calculated by the Beth-Uhlenbeck formula 
using the experimental phase-shifts from [301 ] in compar- 
ison with the excluded volume virial coefficient for dif- 
ferent values of r a . In the triplet state we have a strong 
repulsion between electrons and atoms, hence the hard- 
sphere potential can be expected to give reasonable re- 
sults. Because of the bound state formation, the sin- 
glet state can not be treated within the excluded vol- 
ume approach. Note that in contrast to the excluded 
volume concept and the hard-sphere model, our results 
indeed depend on the temperature. At high tempera- 
tures (T > 50000 K), the second virial coefficient from 
our Beth-Uhlenbeck calculation approaches the excluded 
volume virial coefficient for the atom radius r a = 1.2 as- 
In this sense, the Beth-Uhlenbeck using experimentally 
validated scattering phase-shifts provides a benchmark 
to the semi-empirical excluded volume model. 

Although the excluded volume concept is widely ap- 
plied to take into account the presence of atoms in the 
plasma, this method gives only approximate results. The 
dependence of the atom subsystem on the plasma param- 
eters was included in the confined atom model [1, HH due 
to an atomic radius. These methods cannot cover nu- 
merous effects in the electron-atom interaction, such as 
the spin dependence, scattering phase-shifts, and bound 
states which are included in the Beth-Uhlenbeck formula. 



V. GENERALIZED BETH-UHLENBECK 
APPROACH 

The virial expansion can be extended to higher densi- 
ties if the effects of the medium are taken into account. 
In particular, we outline the consequence of Pauli block- 
ing on the two-particle properties, that is of importance 
when the electrons become degenerate. There are other 
medium effects such as screening, where the effective in- 
teraction potential between the electron and the atom 
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FIG. 8: Second virial coefficient in the triplet channel as func- 
tion of temperature. The Beth-Uhlenbeck result (solid line 
with circles) is compared to the excluded volume model for 
different values of r a 



is replaced by a screened potential. The influence on 
the scattering processes using screened versions of the 
Buckingham and the RDO models has been treated in 
Refs. [ID, [35[ and will not repeated here. A systematic 
approach to screening effects is given within the Green's 
function theory [6|. 

We consider the two-particle effective wave equation 



Pi 



A SE b!) 



4 



+ A SE ( P2 ) 



2mi 2m2 
l±/(pi)±/(pa)] V(piP^PiP2)Wi,P2) 

E(P,T,(i c )tP( PuP2 ), (25) 



Pi,p 2 



where A SE (p) denotes the self-energy shift, and f(p) = 
1/ ^exp [(3(-§^ — m)] T are the Bose and Fermi distri- 
butions. This approach has been applied to charged- 
particle systems as detailed in Ref. [6| for the electron- 
ion system as well as for the electron-hole system. We 
will use a similar approach for the e — a problem. 

The inclusion of self-energy, screening and Pauli block- 
ing effects in the solution of the in-medium Schrodinger 
equation for the electron-ion system leads to non-ideality 
contributions. In particular, the Mott effect is obtained, 
i.e. the dissolution of bound states in the continuum of 
scattering states at increased densities. The contribution 
of the energy shift of atomic levels on the thermodynam- 
ics of the dense hydrogen was considered in Refs [y, [4t| . 
The Pauli shift A Pauli = 32nn e (in Rydberg units) at 
low temperatures and at low densities and the Fock shift 
A Fock = — 207rng, lead to modified behavior at high pres- 
sures. In Ref. [5(|, the effects of Pauli blocking on trans- 
port properties of dense plasma were investigated by solv- 
ing the thermodynamical T-matrix for the electron-ion 



scattering for a separable electron-ion potential. 

A generalized Beth-Uhlenbeck formula has been suc- 
cessfully elaborated for nuclear matter Q . In particular, 
the Mott effect can be included so that the applicabil- 
ity of this approach is extended to the region where a 
quasiparticle description is possible, e.g. in the region of 
strong degeneracy. Analytical expressions are derived for 
a separable potential approach where the in-medium T- 
matrix including Pauli blocking effects can be calculated. 

We study the shift of the binding energy of H~ as 
well as the modification of e — a scattering phase-shifts 
due to Pauli blocking. Our starting point is the effective 
Schrodinger equation for the e — a problem 



L 2m 



l-f(p)\J2v(p,p')^(p') = E il,(p). 
p> 

(26) 

Medium effects are the self-energy (Fock term) A SE (p) 
and the Pauli blocking term [1 — f(p)], that describes the 
occupation of phase space. 

To investigate the Mott effect with respect to the for- 
mation of H~, we investigate the binding energy of the 
e — a system as a function of density, i.e. the differ- 
ence between the bound state energy and the continuum 
edge of scattering states. The self-energy of electrons 
A SE (p) due to the electron-atom interaction shifts both 
the bound state energy as well as the scattering states, 
the net effect on the ionization energy hence being zero. 
The leading term is the Pauli blocking term, that will be 
evaluated in the following. 

We determine the occupation number f(p) in Eq. (126|) 
via the chemical potential /i e according to 



d 3 p 



(2^) 3 ex p [/?(£-- Me)] + 1 



T 



(27) 



Using the parameters of the rank-one separable poten- 
tial given in Tab. I, the binding energies of the electron 
in the negative hydrogen ion have been calculated in de- 
pendence of the total electron density and the tempera- 
ture. The numerical results for the in-medium binding 
energies are given in Fig. [51 We see that the binding en- 
ergy is decreased with increasing electron density. For 
T = 10000 K, at the density exceeding the Mott density 
n Mott _ g g , ^q21 cm -3 ) bound states cannot be formed. 

The influence of the medium on the scattering phase- 
shifts is obtained by solving the T-matrix including the 
Pauli blocking term. Results for different densities are 
presented in Fig. [TU1 At the Mott density we observe a 
jump of the in-medium s-wave phase-shift ryo by n ac- 
cording to the Levinson theorem. 

Using the density dependent phase-shifts and bind- 
ing energies calculated from the in-medium Schrodinger 
equation (|26p , we can calculate the scattering and bound 
parts of the second virial coefficient. 

The results are summarized in Tabs. IIVI and IVl With 
increasing density the bound part of the second virial 
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FIG. 9: H binding energy in dependence of total electron 
number density n e for different temperatures T. 
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FIG. fO: s-wave scattering phase-shifts as function of 
wavenumber k for different electron densities n e at T — 10000 
K. For comparison, the low-density limit (n e = 0) is also 
shown. 



coefficient is decreasing because the binding energy be- 
comes smaller due to the Pauli blocking and screening 
effects. 

It should be mentioned that the account of the self- 
energy A SE (p) would contribute to the chemical poten- 
tial as determined by the normalization condition (|27p . 
If the density of the electrons is replaced by the density 
of quasiparticles that account for the self-energy shift, 
we also have to modify the Beth-Uhlenbeck expression 
for the second virial coefficient as shown in Ref . [8| . The 
Beth-Uhlenbeck expression for the second virial coeffi- 
cient (J7|) used here is consistent with the single particle 
contribution given by free particles as done here. In fu- 
ture investigations, an improved treatment of density ef- 



fect can be performed on the basis of quasiparticles and 
their corresponding screened interactions. 



VI. CONCLUSIONS 

For partially ionized plasmas, the cluster virial expan- 
sion for thermodynamic functions is considered. We fo- 
cus on the contribution due to the electron-atom interac- 
tion. With the help of the Beth-Uhlenbeck formula, the 
second virial coefficient in the electron-atom channel is 
related to phase-shifts and possible bound states in that 
channel. In contrast to former approaches, we give val- 
ues for the second virial coefficient in the e — a channel 
that are not based on any pseudopotential models but 
are directly related to measured data. Depending on the 
accuracy of presently available experimental data, these 
exact results for the second virial coefficient can serve as 
a benchmark to test other more empirical approaches us- 
ing pseudopotentials or related concepts to evaluate the 
thermodynamic properties of partially ionized plasmas. 

From the theoretical point of view, the e — a interac- 
tion amounts to a three-particle problem. At present, the 
most reliable numerical solutions are obtained from vari- 
ational calculations. After comparing these results with 
experimental scattering data, the second virial coefficient 
is presented in the range from 5 xlO 3 to 10 5 K. 

The accurate calculation of the free energy excess due 
to electron-atom interaction is compared with excluded 
volume results that are widely used in the chemical 
model. This semi-empirical treatment contains the hard- 
core radius of the atom as an empirical parameter. Com- 
paring the corresponding virial expansions, it is shown 
that a single parameter choice for the hard-core radius 
cannot reproduce the non-ideal contribution to the ther- 
modynamic functions in a wide region of temperature. 

We also considered different empirical pseudopoten- 
tials that can approximate these microscopic input quan- 
tities. In particular, a rank-two separable potential is in- 
troduced that fits the microscopic data. The advantage 
of a properly chosen pseudopotential is that higher order 
non-ideality terms with respect to the density can be cal- 
culated. For this, the on-shell properties in the two-body 
channel are no longer sufficient. 

Going beyond the second virial coefficient, density ef- 
fects such as self-energy shifts and Pauli blocking have to 
be considered. In particular, we include Pauli blocking 
using the separable potential. In this way, we performe 
calculations for the density dependent second virial coef- 
ficient to cover electron densities up to near solid densities 
n c ~ 10 22 cm~ 3 . 

Our study of the contribution of the electron-atom in- 
teraction is a step to the systematic evaluation of the 
thermodynamics of partially ionized dense plasma where 
artificial parameters such as a hard-core radius are ob- 
solete. The main ingredient, the systematic transition 
from the physical picture to a chemical one, can be ob- 
tained from a quantum statistical approach. The use of 



13 



TABLE IV: Results only for s -wave:the singlet scattering part of the second virial coefficient bJ c e ' slnglet (T, n e ) for e—H interaction 
for different electron number densities and temperatures. The phase-shifts are calculated with the separable model. 



T, K 


n e = 10 18 cm -3 


n e = 10 19 cm" 3 


rie = 10 21 cm -3 


n e = 5 * 10 21 cm -3 


5000 


-0.0505 


-0.0504 


-0.0464 


-0.1012 


6000 


-0.0547 


-0.0546 


-0.0521 


-0.0764 


7000 


-0.0584 


-0.0583 


-0.0571 


-0.0750 


8000 


-0.0617 


-0.0616 


-0.0614 


-0.0770 


9000 


-0.0647 


-0.0646 


-0.0653 


-0.0797 


10000 


-0.0674 


-0.0674 


-0.0688 


-0.0825 


11000 


-0.0700 


-0.0700 


-0.0708 


-0.0851 


12000 


-0.0725 


-0.0725 


-0.0729 


-0.0871 


13000 


-0.0748 


-0.0748 


-0.0755 


-0.0894 


14000 


-0.0770 


-0.0770 


-0.0772 


-0.0917 


15000 


-0.0792 


-0.0792 


-0.0794 


-0.0939 



TABLE V: The bound part of the second virial coefficient t,b° und . sm s lct (y j ne ) f or e _ jj interaction for different electron number 
densities. 



T, K 


n e = 10 17 cm -3 


n e = 10 19 cm" 3 


n e = 10 21 cm -3 


n e = 5 * 10 21 cm" 3 


5000 


4.8698 


4.8298 


2.1955 


0.2505 


6000 


2.9688 


2.9506 


1.6384 


0.2865 


7000 


2.0848 


2.0749 


1.3154 


0.3243 


8000 


1.5993 


1.5932 


1.1080 


0.3529 


9000 


1.3013 


1.2973 


0.9650 


0.3730 


10000 


1.1034 


1.1006 


0.8610 


0.3864 


11000 


0.9640 


0.9620 


0.7824 


0.3948 


12000 


0.8615 


0.8599 


0.7210 


0.3996 


13000 


0.7833 


0.7820 


0.6719 


0.4019 


14000 


0.7219 


0.7209 


0.6318 


0.4023 


15000 


0.6726 


0.6718 


0.5984 


0.4015 



the technique of Green's functions allows for the account 
of higher order many-particle effects. 

The generalization of the cluster-virial expansion, in- 
cluding more bound states as well as excited states, is 
straightforward. 
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Appendix A: Cluster virial expansion and 
thermodynamic functions 

We consider a system of interacting electrons (e, s e = 
1/2), ions (protons i,S{ — 1/2), and atoms (a) in the 
singlet state (s a . singlet = 0) and in the triplet state 
(sa. triplet = !)■ We neglect the hyperfme splitting. 

Within the cluster virial expansion, the grand canon- 
ical partition function, fi(z e , z ii z a, T, V) is a function of 
the fugacities, z c — e^"^^ ') (see Sec. II), the temper- 
ature T, and the system's volume V. We expand up to 
the second order in the fugacities: 

Q(z c ,T,V) = 1+51 z &c{T,V) 

c—e,i,a 

+ J2 z cZ d n c d{T,V) + 0(z 3 c ) . (Al) 

Here, we have introduced the single-particle partition 
functions f2 c (T, V) = (2s c + 1)V/ A^ and the two-particle 
partition functions Q cc i(T, V) that are related to the in- 

1 /2 

teraction, A c = ^2ttH 2 /m c k^T^ is the thermal wave- 
length of species c. The two-particle partition will be 
related below to the second virial coefficient. From the 
partition function, we can directly derive the pressure P 
in the system, 

P(z c , T, V) — In Q(z c , T, V) . (A2) 



Replacing the partition function by Eq. (|Aip and expand- 
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ing again up to second order in z%, we arrive at 



After some lengthy manipulations we obtain the result 



P(z c ,T, V) _ x ^ 2s c + 1 
k B T ~ L> A? Zc 

c <- 



cd 



(A3) 



We have introduced the second virial coefficients b cd = 

A 3 l 

2s c +1 (£l cd — ^Vl c Vl d )/V that are related to the symmetric 

A 3 ~ 

expressions b cd = 2^+l b cd- 

The thermodynamic functions of the partially ionized 
hydrogen plasma are derived from the grand canonical 
potential Eq. (|A1[) . First, we evaluate the number densi- 
ties of each component, 



d lnCl{z c ,T,V) 



8Z C 



V 



d P(z c ,T,V) 
dz c k B T 



T,V 



T,V 



(A4) 



In terms of the second virial coefficient, evaluation of the 
derivatives with respect to the fugacity yields 



2s c + 1 
"AT" 



[z c + 2 £ Zc z d ^j±±b cd + 0(z 3 c )j 



In the first terms of equations (|A5[) , we recognize the 
partial densities of the ideal (non-interacting) system, 
ri Cjid = (2s c + 1) z c /A 3 c . 

Next, we evaluate the entropy density of the partially 
ionized plasma [z c is a function of /i c and T), 

S(z c ,T, V) = f dP(z c ,T,V y 

V \ dT /Zc<v 



e f mn { z c ,T,v) 

dT V V 



z c ,V 



P{z c ,T,V) k B T f dlnQ{z c ,T,V) 



T 



V 



dT 



(A6) 



5P 



- — - k B ^2 n c In z c + k B T ^ 



(2s c + l)z c z d db cd (T) 



dT 



(A7) 

The density of internal energy follows from the relation 

U(z c ,T,V) TS 



V 



V 



H c n c - P 



3„ „ . , m9 v-^ (2s c + l)z c z d db cd 



'-P-J2n c E c + k B T 2 J2- 



2 ' ■ — - A 3 dT 

c cd r \ 

And finally, wc find for the free energy density f(n c , T) 
F/V after eliminating z c with Eq. (|A5[) 



(A8) 



f(n c ,T) 



U TS 

V ~ ~v 

= ^ ^n c - P 



Ah 



= fc fl T^n c [ln^ I -l + /34°)] 

C c 

- k B T^n c n d b cd (A9) 



that coincides with the expression 



Appendix B: T-matrix for separable potential 

We consider the low density limit n e A 3 <g; 1 , where the 
T-matrix equation [5l| for a separable potential reads: 

T(p,p',E) = V(p,p') 

+ E^'P")^ T(p\p,E). (Bl) 

- ^ E 



2m r 



If we consider a rank- two separable potential (|15p . we 
obtain the following expression for the T-matrix: 



f d 3 v" 1 

T(p,p,E) = \ lWl {p) Wl (p') + X 2 w 2 (p)w 2 (p') + \ lWx (jp) / -nr^wxip")— T(p",p',E) 

J (2tt)3 i E 



d 3 p" , „. 1 



+X 2 w 2 (p) I "wzfr") ^ - T(p",p',E) 



2m r 

= ci 1 w 1 {p)wi{p') + c 22 w 2 {p)w 2 (p') + c 12 wi(p)w 2 (p') + c 2 iw 2 (p)wi(p'), (B2) 

I 

with This set of equations can be simplified if we introduce 

f d 3 p" „ 1 the integrals I tJ (E) = J i c Wt ( p '>.(p") : 

Cij = Ai^ + A, y J^E w i(P >—2 -(cijWi(p ) h^- E 

^ ~ Cij = XiSij + Xihx {E)c Xj + XJ l2 (E)c 2j (B4) 

+ c 2j w 2 {p")) (B3) 
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Finally, after some mathematics, we obtain the T- 
matrix in the following form: 

T(P,P',E) = ^±^{(X 1 [l-X 2 h2(E)}w 1 {p)w 1 (p') 

+ A 2 [l- \ 1 h 1 (E)] W2 (p)w 2 (p') 
+ XiX 2 h 2 (E)w 1 (p)w 2 (p') 

+ X 1 X 2 I 21 (E) W2 (p)w 1 (p r )} (B5) 



where 



Det(E) = l-Ai/ii(^) 1-X 2 I 22 (E) 
— \\\ 2 Ii 2 {E)I 2 i{E). 



(B6) 



Using properties of the T-matrix, the binding energy 
Eq can be obtained from the equation Det(Eo) = 0. 



